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Survival probability of a diffusing test particle in a system of coagulating and 

annihilating random walkers. 
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We calculate the survival probability of a diffusing test particle in an environment of diffusing 
particles that undergo coagulation at rate Ac and annihilation at rate Xa- The test particle dies at 
rate A' on coming into contact with the other particles. The survival probability decays algebraically 
with time as t~^ . The exponent S in d < 2 is calculated using the perturbative renormalization 
group formalism as an expansion in e = 2 — d. It is shown to be universal, independent of A', and to 
depend only on S, the ratio of the diffusion constant of test particles to that of the other particles, 
and on the ratio Aa/Ac. In 2-dimensions we calculate the logarithmic corrections to the power law 
decay of the survival probability. Surprisingly, the log corrections are non-universal. The 1-loop 
answer for 9 in one dimension obtained by setting e = 1 is compared with existing exact solutions 
for special values of 5 and Xa/Xc- The analytical results for the logarithmic corrections are verified 
by Monte Carlo simulations. 

PACS numbers: 05.10.Cc, 05.70.Ln, 82.40.Qt 



I. INTRODUCTION 

The calculation of the survival probability of a test 
particle in reaction diffusion systems has been studied in 
different contexts such as sitepersistence Q , walker per- 
sistence problems |^ |^ ^ ll:.!! polydispersity exponents 
in models o f ag gregation [a, S 8, 9] and predator-prey 
models lllj | The approach to these problems have 
mostly been based on studying exactly solvable limiting 
cases or using mean field approximation and its improve- 
ments such as the Smoluchowski approximation Il3j | . 
In recent years, field theoretic methods have proved suc- 
cessful in providing a general framework to understand 
these problems. In particular, the renormalization group 
analysis has been instrumental both in identifying the 
universal persistent features of reaction-diffusion systems 
and in extracting quantitative results about persistence 
exponents which could not be obtained using other meth- 
ods 1^ lli Q [H . In this paper, we apply field 
theoretic methods to the problem of the survival prob- 
ability of diffusing test B-particles in a background of 
diffusing ^-particles undergoing the reactions 
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The above reaction has been studied in the context 
of persistence. In one dimension and when _B-particles 
are stationary, calculating the survival probability of B- 
particles is equivalent to calculating the fraction of spins 
that have not flipped up to time t in the q-state Potts 



model evolving via zero temperature Glauber dynamics, 
where q = -j^ + 2 [J . The more general problem in which 
the i3-particles are mobile with a diffusion constant equal 
to 6 times the diffusion constant of the A-particlcs has 
been studied in Kefs. 0, H E IS- The density of B- 
particles then decays with time as t~^^^'^\ where Q = 
(Ac + Aa)/(Ac + 2Aa). As Q varies from 1/2 to 1, the ratio 
Xc/Xa varies from to cxd. The known results for 6(6, Q) 
are briefly reviewed below. 

When the dimension d is greater than the upper critical 
dimension- two in this case- the decay exponents are ob- 
tained by solving the mean fleld rate equations with ap- 
propriately renormalized lattice-dependent reaction rate. 
In dimensions d < 2, fluctuation effects become impor- 
tant, and 9(6, Q) is no longer given by the rate equations. 
Exact solutions are one way of calculating exponents in 
1-dimension. When 5 = 0, by mapping the calculation 
of the persistence probability to a calculation of empty 
interval probabilities in the A + A A model, it was 
shown that jlj] 
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Attempts to generalize the methods used in Ref. [l| to 
arbitrary S were successful only in determining the values 
of ^|q=o and ^|(5=o Another solvable limit is Q = 
1, when annihilation is absent. In this case, the problem 
reduces to a three particle problem which can be solved 
exactly to yield (l^ 



2C0S-1 [6/(1 + 6)]'' 
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More general ways of understanding the effects of fluc- 
tuation in low dimensional reaction-diffusion systems in- 
clude the Smoluchowski approximation |l2l |l3j which ef- 
fectively reduces to the replacement of the reaction rates 
in the rate equations by diffusion- renormalized values, 
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and the renormalization group formalism. The exponent 
6(6, 1/2) was calculated using the Smoluchowski approx- 
imation in Ref. The advantage of Smoluchowski 
approximation is its computational simplicity. However, 
it is not clear how to improve this approximation in a 
systematic manner. Also, it was shown in Refs. 
that while this approach gives an answer close to the ac- 
tual one for Q — 1/2, it becomes increasingly worse as 
Q nears 1. The field theoretic approach using the renor- 
malization group formalism currently provides the only 
systematic way of calculating the decay exponents below 
the critical dimension. The exponent 0((5, 1/2) was calcu- 
lated using field theoretic methods in Ref. However, 
the renormalization group scheme used was complicated 
and did not capture the right logarithmic corrections (see 
Sec lIVIBl for a more detailed discussion). 6(1, Q) was cal- 
culated as an expansion in (2 — d) in in the context 
of domain wall persistence in the Potts model. 

In this paper, we extend the formalism developed in 
Refs. 6, 9] to calculate 9{5, Q) for arbitrary 6 and Q to 
order e, where e — 2 — d. In particular we show that 
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where 



f{6) = l-25 + 2(51n 



1 



du 



ln(l - u) 
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The function f{5) increases from (1 — 7r^/4) to as J 
increases from to oo. In 2-dimensions, we calculate 
logarithmic corrections to the power law decay and show 
that 

(6) - ln(^)"/^ (6) 

where (&) is the mean density of B-particles and 



a = Q(l + (5) 



2, + Q{l + 5)f {5) + 2\n- 



A' (1 + ^)(A„ + A,) 



, (7) 



with the function / as defined in Eq. A surprising 
feature of Eq. Q is its non-universality for finite reac- 
tion rates, A, A' < oo. In this case a explicitly depends 
on both reaction rates. This is contrary to the usual be- 
lief that below the upper critical dimension the kinetics 
is diffusion limited and hence one may set reaction rates 
to infinite. Most exact solutions make use of this simpli- 
fying assumption. The above result serves as a counter 
example. 

The rest of the paper is organized as follows. In Sec. ITU 
the model is defined. In Sec. IIIII the rate equation ap- 
proach is compared with the Smoluchowski approxima- 
tion. The survival probability is calculated to one loop 



precision. In Sec. IIVI the renormalization group analy- 
sis of the problem is carried out and Eqs. 10), and 
(|7|l are derived. The one-loop answer for 9 is compared 
with the result of Smoluchowski approximation and also 
with known exact results in 1-dimension for special val- 
ues of 5 and Q. We also compare the analytical results 
with the results from numerical simulations. First, the 
predictions for the logarithmic corrections to the power 
law decay are confirmed numerically in the limit of in- 
stantaneous reactions. Second, the non-universality of 
logarithmic corrections for finite reaction rates is veri- 
fied. Finally, we end with a summary and conclusions in 
Sec.lVl 



II. THE MODEL. 

Consider a d-dimensional hyper-cubic lattice whose 
sites may be occupied by j4-particles and i3-particles. 
Multiple occupancy of a site is allowed. Given a con- 
figuration of particles, the system evolves in time as fol- 
lows, (i) At rate D/{2d), an ^-particle hops to a nearest 
neighbor site, (ii) At rate SD/{2d), a B-particle hops to 
a nearest neighbor site, (iii) At rate Aq, two A-particles 
at the same site annihilate each other, (iv) At rate Ac, 
two A-particles at the same site coagulate together, thus 
reducing the number of A particles by one. (v) At rate 
A', a i3-particle is absorbed by an j4-particle at the same 
site. The initial number of A- {B-) particles at the lattice 
sites are chosen independently from a Poisson distribu- 
tion with intensity oq (6o)- 

The action corresponding to the continuous limit of 
the model can be derived from the master equation using 
Doi's formalism [Hill 113. We skip the derivation and 
give the final result. The action is 

S ^ j dt j d'^x {a{dt - V2)a + h{dt - 5V^)b + ^aa^ 

+ ^a^a^ + A'&a6 + X'abab - {dao + bbo)S{t)^ (8) 

where the a and b are complex fields corresponding to A- 
and _B-particles, the diffusion constant D has been set 
equal to 1, and 



Q = A±^, 

Ac -f 2Aa 
A — Aa + Ac. 



(9) 
(10) 



The knowledge of all the correlation functions of the fields 
a, b allows one to construct all the correlation functions 
of local densities of A and B particles In particular, 
the mean density of A and B particles at {x, t) is equal 
to {a{x,t)) and {b{x,t)) respectively, where (. . .) denotes 
functional average with respect to the functional measure 
Eq. ©. 

The action can be brought into a more convenient form 
by rescaling the fields as follows: {a,b) Q^^{a,b), 
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(ao,6o) ^ Qiao, bo), and (A, A') ^ 2Q(A,A'). Then 
S = J dtj d'^x (^d{dt -\7^)a + b{dt - S\/'^)b + \aa^ 
+ Xa^c? + 2X'Qbab + 2X'abab - [dao + 66n)5(t) 



The Feynman diagrams corresponding to the action in 
Eq. are shown in Fig. ^ 

We are interested in the mean density of _B-particles 
in the hmit of large time, as the survival probability is 
proportional to mean density. The evolution of mean 
density of yl-particles (a) is independent of the statistics 
of i?-particles and decays at large times t as 0| 



(a) 



d<2, 



<-Mn(t), d = 2, 
t-\ d>2. 



(12) 



III. 



THE COMPUTATION OF PERSISTENCE 
EXPONENT USING MEAN FIELD AND 
SMOLUCHOWSKI APPROXIMATIONS. 



The perturbative expansion of (6) in powers of A can 
be constructed using the Feynman diagrams shown in 
Fig. n^]. Diagrammatically, (a) ((6)) is the sum of 
all Feynman diagrams with one outgoing a-(6-)line re- 
spectively. Clearly, there is an infinite number of dia- 
grams contributing to (a) and (6). These diagrams can 
be grouped together according to the number of loops 
that they contain, thus giving rise to the loop expansion. 
Let e = 2 ~ d. A simple combinatorial argument shows 
that the contribution from a diagram with n loops is pro- 
portional to where g{t) = Xt'^"^ When e < 0, 
the main contribution to (a) and (b) comes from prop- 
erly renormalized tree level diagrams (diagrams without 
loops). When e > 0, the loop expansion fails since for 
large times g{t) is no longer a small perturbation param- 
eter. We therefore conclude that 2 is the upper critical 
dimension. For d < 2 we will use the formalism of pertur- 
bative renormalization group to convert the loop expan- 
sion into an e-expansion and calculate scaling exponents 
as a series in e. 



A. Tree level diagrams 

In d < 2 and small times, tree diagrams give the main 
contribution to the survival probability. Let (a)mf and 
{b)mf be mean field densities given by the sum of contri- 
butions coming from tree diagrams with a single outgoing 
a- line and 6-line respectively. We denote {a)„^f and (&)„if 
by thick solid lines and thick dashed lines respectively. 
The integral equations satisfied by (a),„f and (&)rnf are 
presented in diagrammatic form in Fig. |2Ia) andj^^b). 
After differentiating with respect to time, they can be 
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FIG. 1: Propagators and vertices of the theory. 



(a) 



-X = 



(b) X 



2\'Q ^ 



(c) 




(d) 



FIG. 2: Diagrammatic form of mean field equations for (a) 
mean particle density (b) , (b) mean density of _B-particles (b) , 
(c) Gl^, and (d) Gil 



written in analytic form as 

dt{a) = -X(a)\ 
dt{b) = -2QX'{b){a), 



(13) 
(14) 



in which one can easily recognize the rate equations of the 
model. Thus, the identification of tree-level truncation 
with mean field approximation is justified. 

Equations (|13|l and H14() are easily solved yielding 

ao 



(a(0>mf 



1 + Aflgt ' 

bo 



(1 + XaotyQ^-'l^ ' 



(15) 
(16) 



where oq and fog are the initial densities of A- and B- 
particles respectively. Thus, 



e{5,Q) = 2Q^, d>2. 
A 



(17) 



The result is explicitly dependent on A', A while being in- 
dependent of S and describes the reaction-limited regime 
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of the problem. It should be mentioned here that the 
above result is valid only in the limit when the reac- 
tion rates are the smallest parameters in the problem, 
i.e. A, A' <ti Iq '^, where Iq is the lattice spacing. In the 
other limit when the lattice spacing is the smallest pa- 
rameter in the problem, the exponents gets modified to 

m 

e{6,Q) ^Q{l + 6), 1^-^<$:X,\', d>2. (18) 

In order to estimate the validity of the mean field ap- 
proximation in d < 2, the one-loop corrections to the 
mean field answer have to be evaluated. In calculating 
loop corrections to Eqs. IjlSI) and (|16|l . we are faced with 
the problem of summing over infinite number of diagrams 
containing a given number of loops. This problem can be 
simplified by introducing mean field propagators which 
are sums of all tree diagrams with one incoming line and 
one outgoing line. Expressed in terms of these mean field 
propagators, there are only a finite number of diagrams 
with a fixed number of loops. 

Let and be mean field propagators. The 

integral equations satisfied by them are presented in dia- 
grammatic form in Figs.[2Ic) and|2Id). The solutions to 
these equations are 



GNN(2|1) = 



mf J 



(a(ii))mf 
(a(i2))mf 



(a(ii)>r 



Gi(2|l), 

2QX'/X 



(19) 

Gs{2\l), (20) 



where 1 = {xi,ti), 2 = {x2,t2) and Gd is the Green's 
function of the linear diffusion equation with diffusion 
constant D. 



B. Smoluchowski approximation 

Before presenting the renormalization group calcula- 
tion of 9(6, Q), we briefly discuss a method commonly 
used to study fluctuation effects in reaction-diffusion sys- 
tems, namely the Smoluchowski approximation. The es- 
sential idea of the Smoluchowski approach is to relate 
the reaction rates A and A' to the diffusion rates. One 
assumes that particles react instantaneously when they 
come within a fixed radius of each other (see 0, 0| for 
a more detailed discussion). Knowing the first return 
probabilities of random walks, one obtains for c? < 2, 



A - const X t'^/^"^ 



A' ~ const X 



1 + .5 



d/2 



In d = 2 additional log corrections appear and 
const 

A 



A' 



const X (1 + (5) 
21n[(i±^)t] 



(21) 
(22) 

(23) 
(24) 



Replacing A and A' in Eqs. ifT^ and f(T^ by the effec- 
tive reaction reaction rates, and solving for (6), we obtain 



d< 2, 
d= 2. 



(25) 

where {b)s denotes the mean density of i3-particles ob- 
tained from Smoluchowski approximation. In particular, 
the Smoluchowski theory's prediction for 9 is 



9s{S,Q) = dQ[ , d<2. 



(26) 



This answer for 9 depends on 6, Q and space dimen- 
sionality d. It does not however depend on A and A'. 
Thus, unlike the mean field answer Eq. (|17|l . it has the 
correct universality properties for a quantity describing 
reaction-diffusion systems in the diffusion-limited regime. 
However, the Smoluchowski result differs considerably 
from the correct result when Q nears 1. For example, 
in one dimension 9s{l, 1) = 1.0 while 9{l, 1) = 1.5 (see 
Eq. Q)- For more comparisons, see Sec. rrvTI It is not 
clear how one can improve the Smoluchowski approxi- 
mation. The renormalization group formalism, though 
more involved, provides a systematic way to go beyond 
the Smoluchowski approximation. 



C. One loop diagrams 

The rate equation results do not depend on the dif- 
fusion coefficients of the particles or the dimensionality 
of the ambient space. These parameters appear in the 
one-loop corrections to the tree level answers. Using the 
mean field propagators and densities, it is easy to classify 
all the one loop diagrams contributing to (&). These are 
shown in Fig. El Skipping the computations we present 
the final answers for contributions corresponding to each 
of the Feynman diagrams in the limit ao oo: 



(a) 
(c) 



32gA'26ot'/^ 



A(aoAt)2Q^'/A[47r(l + d)]'^/^e^ie + 2) ' 
8Q2A'2 6o(i + S)t'/^ 



(27) 



A(aoAt)2QA'/A[47r(i + J)]<i/2 
-256QA'&oi'/^ 



(aoAi)2Q^'/A(8^)d/2e2(g 2)^{e + 4) ' 



[/(<5) + 0(6)](28) 



(29) 



where (a), (6) and (c) refer to the contributions from 
diagrams in Figs. |2fa), |3Ib) and|3{c) respectively, and 



f{5) = 1-2(5 + 2(5 In 



1 + 5 



+ {1-6') 



du 



ln(l 



(30) 

Adding these one-loop contributions to the mean field 
answer Eq. (|16|l . we obtain in the limit oq — s- oo. 



m) 



A 



l2Qy/\ 



1 



4A' 



1 



8gA'i^/2 

(47r)^/2e \ A(l + 6^/^ e(e + 2) 



5 



(a) 




(b) 




FIG. 3; One loop corrections to the mean field result for (b). 



32 



QA'(l + 5)^/2/(<5) 



2d/2g(£ + 2)2(e + 4) 2A 
2- and higher loop corrections, 



(31) 



where A = bo/iaoXf^^' . We see that if A A', then 
the mean field answer Eq. (|17|l is correct in d < 2 only if 
QAt'^/2 <ti 1. Clearly, this condition breaks down in the 
limit of large times in d < 2. 



IV. PERTURBATIVE COMPUTATION OF 
eiS,Q) NEAR d = 2 USING RG METHOD. 

In this section, we calculate the large time behavior of 
(6). The loop expansion for {b{t)) fails at large times in 
d < 2. To extract the large time behavior of {b{t)) in 
d < 2 we will use the formalism of perturbative renor- 
malization group. 

The renormalization group formalism used in Refs. 0, 
|9[ for the case 6=1 was based on the Callan-Symanzik 
equations for the mean density of _B-particles. There were 
two relevant couplings for the theory: the reaction rate 
A and the initial density Bq. The anomalous dimension 
of {b{t)) was attributed to the renormalization of bo. 

This approach turns out to be very complicated when 
6^1. This is due to the explicit dependence of the 
classical scaling dimension of {b{t)) on A and A'. Further 
complications arise due to non-commutativity of e — > 
and oo — ^ oo limits, which leads to an apparent order- 
l/e^ singularity in the one- loop correction to {b{t)) [see 
Eq. 

These problems are circumvented by analyzing the 
large time asymptotic behavior of the logarithmic deriva- 
tive of {b{t)), rather than {b{t)) itself. It follows from 
Eq. (EH) that 



tdtMm))] = n(i), 



(32) 



where 



n(t) = -2g- 

2-'^/232 



A' AQXt'/^ r 4A' 



(47r)''/2 \X{1 + S)'^/^ e{e + 2) 

, QX'ii + sy/^ns) 



e(e + 2)2(e + 4) 



2A 



O(A^). (33) 



The large time asymptotic behavior of 11 (i) can be ob- 
tained by solving the Callan-Symanzik equation with ini- 
tial conditions given by the right-hand side of Eg. 
regularized at some reference time to (see Ref. |22 for 
a review). The coefficients of Callan-Symanzik equation 
are determined by the law of renormalization of all the 
relevant couplings of the theory. Power counting analo- 
gous to that carried out in Ref. 9] shows that there are 
only two relevant couplings of the theory which deter- 
mine large time behavior of 11 (t) in d < 2: the reaction 
rates A and A'. We mention here that ll{t) is simply 
related to the polarization operator used in Ref. 0. 



Let go = Xt\ 



9o 



Xtn^ be the dimensionless reac- 



tion rates. We choose in such a way that go, g'o 1- 
The law of renormalization of reaction rates has been 
worked out in Ref. 0, 0|. The renormalized reaction 
rates gj^ and g'j^ are related to go and as follows: 



9B. = 
9r = 



90 



1 + 90/9* ' 
9'o 

l+9o/9'.' 



(34) 
(35) 



Here and are the nontrivial fixed points of the renor- 
malization group flow in the space of effective coupling 
constants, and are given by 



9* 
9'* 



2r(e/2)' 
[4^(l + ^)]''/2 
2r(e/2) ' 



(36) 
(37) 



where T{x) is the Euler Gamma function. The renor- 
malization of both coupling constants is due to the same 
physical effect - the recurrence of random walks in d < 2. 

n(io) expressed in terms of g^ and has the following 
form: 



n(io) = 



Q9'r 



9R 
Q9'Hf{5) 
29R 



5fl(7"l) 5-27 



0{e) 



0{g' 



(38) 



where 7 is the Euler constant. n(io) regarded as a func- 
tion of gn and g'j^ is non-singular at e = 0. As a result this 
expression is valid for d < 2. The lack of io-dependence 
of n(i) for i > to is expressed by the following renormal- 
ization group (Callan-Symanzik) equation: 



■ d pjgR) d (3{g'j,) d , 
. 9i 2 dgn^ 2 dg'j,' 



where /3(gi?) = -2todgBjdto and /3(5^) 
are the beta functions. 



n(t) = 0, (39) 
-2todg'Jdto 



K9r) 



9r{9r - 



(40) 
(41) 
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We will now solve Eq. (|39|) with the initial condition given 
by Eq. 138() to obtain the large time asymptotic behavior 
of n. We then extract the large time asymptotic behavior 
of {bit)) by solving Eq. 



A. Survival probability in d < 2 

At large times, the solutions of Callan-Symanzik equa- 
tion (|39|l are governed by the stable fixed points of the 
beta functions. In d < 2, these are gn — g<, and g'j^ — g'^. 
It then follows that 



n(t) = -2Q{\ + 5) + eQ{l + 5) 
1 + 5 



In 



1 



6 3 
-+2 



Q- 



-f{s) 



(42) 



Substituting Eq. into Eq. and solving for 

{b{t)), we obtain the 0{e) result for 9: 

where the function f{d) is as in Eq. H3U|I . 

We now compare the 1-dimensional result obtained 
by putting e = 1 in Eq. (|43|l with exact results in 1- 
dimension for special values of 6 and Q. The exact result 
for 9{0, Q) in 1-dimension is given in Eq. |5J, while that 
for 9{S, 1) is given in Eq. ©. Figs 0] and |S1 show the re- 
sults for S — and Q — 1 respectively. The e-expansion 
result is seen to compare very well with the exact re- 
sult. On the other hand, the Smoluchowski results fail 
for Q larger than 1/2 when 6 = 0, and for all 6 when 
(5 = 1. It should be noted that when S becomes large, the 
e-expansion should fail and Smoluchowski result should 
give a much better approximation. 



0.8- 



0.6- 



^0.4- 
0.2 - 



0.2 



. 4 



. 6 



. 8 



Q 



FIG. 4: The one-loop answer (Eq. 1431 ') is compared with 
exact result in 1-dimension when 5 = (Eq. 0)- The solid 
line corresponds to the exact answer, dashed line to one-loop, 
and dot dashed line to the Smoluchowski result (Eq. 11251 '). 




FIG. 5: The one-loop answer (Eq. (HHJ) is compared with 
exact result in 1-dimension when Q = 1 (Eq. Q). The solid 
line corresponds to the exact answer, dashed line to one-loop, 
and dot dashed line to the Smoluchowski result (Eq. 1)25^'). 



B. Survival probability in d = 2. 

The upper critical dimension of our model is two. The 
non-trivial fixed points of the /9-functions Eqs. H40|l and 
(|41|) vanish at d = 2. It is then expected that the rate 
equation results give the correct large time behavior of 
mean densities, perhaps modulo logarithmic corrections. 
This turns out to be incorrect. The complication comes 
from the fact that 9 predicted by the mean field theory 
(see Eq. (|17|) ') is non- universal and depends on the ratio 
of coupling constants g^ and g'j^. Each of these couplings 
is marginally relevant in two dimensions and flows under 
RG transformations to as [ln(t)]~^. Their ratio flows 
to a finite universal value which determines the algebraic 
decay of the survival probability. However, the deviation 
from this universal value vanishes with time as C/ln(<), 
where C is a non-universal constant. This results in non- 
universal logarithmic corrections to the universal power 



law decay of survival probability. It is worth mentioning 
that in ci < 2, similar considerations lead to a conclu- 
sion that the amplitude of corrections to scaling is non- 
universal. 

We need to solve the Callan-Symanzik equation H39|) 
with coefficients calculated at d = 2: 



/3(3)U=2 
/3(5')U=2 
Then Eq. reduces to 

9l d , 



El 

2tt' 
9''- 



n{l + 6)' 



d 



'2 

5fl 



d 



(44) 
(45) 

n(i) = 0, (46) 



AndgR 2^{l + 5)dg'^j 
which has to be solved with the initial condition given by 
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FIG. 6: The variation of the mean density of _B-particles 
in two dimensions with time is shown. The simulations were 
done on a 3200 x 3200 lattice for Q = 0.5 and 5 = 0. The data 
has been averaged 1000 times. In the inset, the variation of 
{b)y/t with time is shown. The power of the logarithm in the 
best fit is 0.23 ± 0.03. 
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FIG. 7: The variation of the mean density of B-particles 
in two dimensions with time is shown. The simulations were 
done on a 3200 x 3200 lattice for Q = 1.0 and 5 = 0. The data 
has been averaged 1000 times. In the inset, the variation of 
{b)t with time is shown. The power of the logarithm in the 
best fit is 0.08 ± 0.04. 



Eq. (|2HI) at t = to, e ^ 0. The solution is 



n{t) = -Q{i + s) 



+ tt{1 + 6) 



9'r 



2Q{l + 6) 

Ht/to) 
2 

(1 + S)gR 



Q 



m 



o 



In (^) 



(47) 



The non-universal term in Eq. H47(l is proportional to 
— 2/[(l + S)gii]. It is convenient to express this 
amplitude in terms of bare couplings. In two dimensions, 



1 2 1 



9r (1 + S)gR g'o 



ln(l±i) 



(l + (5)(7o 27r(l + 5)' 



(48) 



In d < 2, Eq. H48I) has to be modified by the omitting the 
logarithmic term on the right hand side. 

Solving Eq. with Eq. substituted for the right 
hand side and taking H48|l into account, one finds that 



(b{t)) = const X 



{gRMt/t,)r^^ 
(g«i)Q(i+^) 



l + 0( 



Ht/to) 



(49) 



where 



a = Q{l + S) 



3 + Q(l 



-47rQ(l + <5)2 



2 



21n- 



il + S)\ 



(50) 



Thus, in two dimensions, the power law exponent is uni- 
versal and independent of A and A'. However, the loga- 
rithmic corrections are non-universal and depend on mi- 
croscopic reaction rates A' = g'o and A = go- However, 
most simulations are done in the limit when the reactions 
are instantaneous. In this limit, the non-universal term 
in Eq. (|50|l is zero. 



A, = ~, X-- 
X = 8.0, X -- 
X = 8.0, X -- 



16.0 
8.0 
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FIG. 8: The variation of the mean density of B-particles in 
two dimensions with time is shown. The simulations were 
done on a 900 x 900 lattice for Q — 0.5 and 5 = 0. The data 
has been averaged 1000 times. 



The log corrections in Eq. (|50|) arc different from the 
log corrections calculated for the Q — 1/2 case in Ref. i. 
This discrepancy is due to the fact that only renormalized 
tree level computations were done in 0, while to obtain 
the correct logarithmic dependence, one-loop corrections 
have to be taken into account. 

We also mention that if one were to ignore the contri- 
bution from one- loop diagrams, then the log corrections 
would be identical with the log corrections obtained from 
the Smoluchowski approximation [see Eq. (|25() ]. and will 
be different from the log corrections obtained from renor- 
malized tree-level as in Ref. 01 • 

We now study logarithmic corrections numerically. 
First, consider the case when the microscopic reactions 
are instantaneous, i.e., A = A' = oo. In this limit, the 
non-universal term in Eq. H5U|) is equal to zero. The 
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Monte Carlo simulations were done for this case on a 
two dimensional lattice of size 3200 x 3200 with peri- 
odic boundary conditions. As the reactions are instanta- 
neous, the maximum number of particles at a site is one. 
The simulations were done for 6 = 0, i.e., immobile B- 
particles. The results for Q = 0.5 and Q = 1.0 are shown 
in Figs.iniandHrespectively. For Q = 0.5, a ^ 0.23±0.03 
which compares well with the theoretical value of ap- 
proximately 0.22 in Eq. jSUJ. Note that renormalized 
tree-level answer is a = 0.5 2J. For Q = 1.0, the nu- 
merical value is .08 ± 0.04, which compares well with the 
theoretical value of approximately 0.07. This answer also 
deviates significantly from renormalized tree-level value 
of 1.0. 

Second, we study the logarithmic corrections for finite 
reaction rates to test the non-universal term in Eq. (|50|l . 
In this case the lattice size is 900 x 900. The results for 
d = 0, Q — 0.5 and different reaction rates are shown 
in Fig. |S1 If A = 16, A' = 8, then the non-universal 
term in Eq. (|50|) is zero and a ~ 0.22 as in the case 
of infinite reaction rates. This is consistent with dashed 
line in Fig. |S| being parallel to the solid line. If however, 
A = 8, A' = 8, then a « —0.57. The slope of the bottom 
line in Fig. |Slis clearly negative. 

V. SUMMARY AND CONCLUSIONS 

In summary, we calculated the large time behavior of 
the survival probability of a test particle in a system of 
coagulating and annihilating random walkers in d < 2. In 
one dimension, this generalizes the site persistence prob- 
lem in the q-state Potts model evolving via zero temper- 
ature Glauber dynamics. The survival probability was 
shown to decay as a power law. In d < 2, the exponent 
9 characterizing this power law was shown to be univer- 



sal, in the sense that it depends only on d and Q and was 
independent of A, A'. The renormalization group formal- 
ism provided a systematic way of calculating the survival 
probability for the entire parameter space. 

In two dimensions, we computed the logarithmic cor- 
rections to the power law decay. It was shown that to 
compute the correct logarithmic factors, one had to in- 
clude contributions from one-loop diagrams and not just 
the tree level diagrams as was done in earlier work. The 
behavior of the survival probability in two dimensions 
is surprising. First, the power law decay is universal 
and thus does not coincide with the rate equation re- 
sult, even though d = 2 is the upper critical dimension. 
Second, the logarithmic corrections to the power law are 
non-universal and depend on the reaction rates. This 
is contrary to the general expectation that kinetics of 
reaction-diffusion systems are diffusion limited below the 
upper critical dimension. Both the universality of the 
power law and the non-universality of log-corrections in 
two dimensions can be traced to the fact that the rate 
equation exponent is given by the ratio of microscopic 
rates, which are both marginally relevant in two dimen- 
sions. 

From the renormalization group point of view, by 
studying the logarithmic derivative of the mean density of 
-B-particles, we have considerably simplified the schemes 
used in Refs. 0jlEPl- While the exponents are only cal- 
culated up to first order in e, the renormalization group 
method remains the only systematic way of computing 
the exponents when an exact solution is not available. 
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